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I. Maxwell Stress Tensor 



A. Notation 



F x= V ^x- T x = T xx T *+T xy ! y +T xz L 



F v = V *V \ = T yx T * +T yy T Y+ T yz T ; 



F z = V-t z , t z = T zx Tx+T zy T y +T zz L 



T = 



"^"xx "^"xy xz 
yx ^yt yz 
~^~zx ~^"zy ~^"zz 



fx = |F x dV = jv.^ x dV=f x x .nda = j [T xx n x + T xy n y +T xz n z ]da 



T x * n = T xx n x + T xy n y + T x Z n z = T x n n n 



\ * n = T yx n x + T yy n y + T y Z n z = T yn n n 



\ ' n = T zx n x + T zy n y + T zz n z = T zn n n 



f, = f V • T,dV = j> t, • n dV = j> T^dS = J F,dV 



F. = V • t = — T ix + — T iv + — T i7 

1 ex lx ay Iy dz lz 



Sx H IJ 



B. EQS Stress Tensor 



F = pE--E-EVs + v(-E-E — p 

f 2 1,2 a P F 



= V.(8E)E-i(E.E)vs + vf^E.E 



58 

a P 
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_ 9 ( eE j) c i c ,. ds d fi_ _ as 



VxE=0^> 



2 * * ex, ex, i^2 

9E, _ BE i 



dx i ex; 



F = — (eE,E,) - eE, - ^-E K E K — - t K t, — p 

9x ; v 1 ' y 3 dx j 2 5X; 5X; 1,2 c3p 



8E, 1 



58 9 f 1 98 



5 ( C ^ \ ^ lr r OS 9 f 1 ._ _ 98 ^ 
F i = SE i E 1 - SE 1 ~ ~ ^ E K E K + ^ E K E K P 

1 9x ( v 1)7 , 9Xj 2 5x 5x ^2 9p J 



d ( 1 _ r 

8 — E i E , 

9x, I 2 J J 



F, = -^-(eE.E.) 

1 axj v 1 >' 



d 

9x ; 



1 ,_ 1 98 _ ,_ 

2-s e k e k ~ 2 P 9p K K 



— = 5 — 

9Xj ij axi 



8 -J° ^ j 



Kronecker Delta 



F, = 



9x f 



eE i E i -l8 n E K E K [e-p — 
1 J 2 J ^ 9p 



— — — 
9x f l] 



C. MQS Stress Tensor 



F = J f x B - -H • H Vu + V 
2 



ip^H.H 

2 K 9p 



= (VxH)x(nHj-^H-H Vp + V^|p|^H-H 

(vxh)xH = (R.v)h-^v(R.h) 



F = n 



(h-v)h--v(h-h) 



-H'HVu + V f-p^-H • H 
2 [ 2 V dp 
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F, = n 



H,— H, 

J dx j 1 



l_d_ 
2 8x, 



(H K H K ) 



1 M M d V 

— n^n,. 



8 ( 1 Su. 



sx, ax : 1^2 ap 



H„H|, 



ax 



-(nH,H J )-H,-|-(nH j ) 



u. 3 



d\x S [ 1 8(i 



2 - x ' K ' ' K 2 ' ' K ' ' K ax : 1 ax : (2 p ap H|<Hk 




F, = 



ax 



— — — ( — uhLH,, - p — HLH,, 

aXj 2 ap 



ax jL 



- -8 ij H K H K 



H - p 



ap" 



-*-T.. 

a Xj IJ 



Ty - (iH^ - -8 ij H K H K 



H - p 



ap" 



II. Air-Gap Magnetic Machines 




staior 




Fi S , 



4 2 L Typical "air-gap" configurations in which a force or torque on a rigid rotor r 
'from spatially periodic sources interacting with spatially periodic excitations on a r 
"stator." Because of the periodicity, the forte or torque car. be represented n_n terms 
electric or magnetic stress acting at the air-gap surfaces (a) planar geometry or 

oped model; (b) planar or cylindrical beam; (c) cylindrical rotor. 



esults 
igid 

of the 
devel- 



Courtesy of MIT Press. Used with permission. 
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A. Generalized Description 



d -- 



o m 



/ / (stater)' y / / / 



K s = Re[K s e-jkz ]iy 



/ 




K r = Re[K r e-i^Jiy 



/ (rotor) JJr* 1 ^ ! / 

/ / 



Depth w 



B.C.: H z (x=A) =K S 
H z (X=0) =-K r 



2u/k 



2ji/k 



f z = J T zx n x dzdy = w j u H z H x dz = w j u H z H x dz 

^ 

force on a wavelength 



i (z, t) = Re [Ae jkz ] , b (z, t) = Re [Be jkz ] 



2 it A 



2ti 



j a(z,t)b(z,t)dz = -Re 



AB 



= -Re 



A B 



f, = 



27iw n 



k 2 



Re 



/^7tW(X 



Re 



H Z H X * 



= n k 



-coth kd 



sinh kd 



sinh kd 
coth kd 
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1 jk jk 



sinh kd 



+ 1 coth kd 



!v - ^— coth kd 



jksinhkd jk 



Re 



-K r H x n 



= -Re 



K K 

r 

sinh kd 



■ + K K coth kd 



= Re 



-|i jK r K s /sinh kd 



f *W n„ Re 

k sinhkd 



(force on each wavelength) 



B. Synchronous Interaction 




Fig. 4.3.1. Rotor and stator reference 
frames z' and z. 



Courtesy of MIT Press. Used with permission. 
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K s = K* sin[co s t-kz]= Re[-jK=e j(Mst " kz) ] 
K r = K r sin[co r t-k(z'-8)]; z' = z-Ut 

= K r sin[(co r +kU)t-k(z-S)] 

= Re[-jK r e j(l0r+kU)t e jk5 ] 

K = -jK= e Kt 
K = -jK r e jk5 e jK+kU)t 

For time average force => ro s = ro r + k U (synchronous condition) 
Usually ra r = => a> s = kU 



(f ) = .*JLK« K^sink5 
w k sinhkd 



Table 4.3.1. Basic configurations illustrating classes of electromechanical 
interactions and devices. MQS and EQS systems respectively in 
left and right columns. 




Courtesy of MIT Press. Used with permission. 
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III. Electrostatic Machine 



/ 

d - 


I* 




O s = Re[V s e-J kz ] 








C 


& > 


- — 2rc/k — ^ 


O r = Re[V r e-i 10 ] 


1 z 


+ + + 


Depth w 



f = 



w2ji 



2n/k 



dz = 



2jiw 



2jt/k 



x=0 



J B E Z E X 



x=0 



E z = jkV r 



z X k 



s E r E r 

z x 



TtW 



Re 



^(-jkv;)E x 



D x 
D x 



s k 



-coth kd 
1 

sinh kd 



sinh kd 
coth kd 



Vs 
Vr 



£ E x — s k 



-V 



sinh kd 



+ V coth kd 



Re 



"jks V r E x 



= Re 



-jk 2 s V r 



-V 

+ V coth kd 

sinh kd 



= Re 



+jk 2 8 V s v7sinhkd 
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f - " w i<^o Pc 

z X sinhkd 



iV V 

■> s r 



V s = V S cos (oo s t - kz) 

V r = -V r cos (ro r t - k (z' - 5)) ; z' = z - Ut 



V =-V r e 
V s = V r e Kt 



r„j(%+kU)t a jk8 



jiwks 
sinh kd 



Re 



[-J Vc K 



g-jkS e j(m s -co r -kU)t 



ro s = co r + kU 



sinh kd 



IV. Derivation of the Korteweg-Helmholtz Force Density for Incompressible Media from 
the Quasistatic Poynting's Theorem 

A. Poynting's Theorem 

3B 



VxE = 



at 



VxH = Jf + 



3D 
~dt 



V-D = Pf 



V«B = 



V.(ExH) =H.(VxE)-E.(VxH) 



5t 5t f 
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B. Power In Quasistatic Electric Circuits 



Vn-i 


\ 


h 
h 


,v 2 

/l2 Vi 










h 


\s 




*S= E x H 


The circuit power into an N terminal network 
Sk-i V A equals the electromagnetic power flow 
into the surface surrounding the network, 
-cf s ExH.dS 



Far away from the circuit elements 

VxE = 0=>E = -V<S 

VxH = J f ^V.Jf = 



P in = -^ExH).da 

s 




= -J f .v<D = -V.(j f <X>) 
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R n =-fv.(j f <D)dV 

V 

= -j> Jf <6 «da 



= -2>k #J f -da 

k=l s 



-L 



= Iv k i k 



C. Electroquasistatics (EQS) 
Ohmic Media: J f ' = ctE ' = J f - p f v => J f = ctE + p f v 

D = e(x, y,z)E 



J"v.(ExH)dV = | Ex H- da = -£V k I k = -j" E * J:( £ ( x 'y' z ) E ) dv 



-j*E.(aE + Pf v)dV 

v 

^ y ^=l E -^^'U^ Z ^) dV + jHE| 2 dV + J Pf E.vdV 



I e(x,y,z) St 
J v 2 s(x,y,z 



s 2 (x, y,z) E 



J 



1 5 



J v 2s(x,y,z) St 



)stL 



s 2 (x,y,z)|E| 2 dV = J 



St 



dV + Jcr|E| 2 dV + J Pf E-vdV 

v 

s/(x, y,z)ir 



-J 



s 2 (x,y,z)|E| d 



St 



dV 



' 1 A 
s(x, y,z) 



dV 



-J 



■e(x,y,z)|E|" 



dV + jU_|(s(x,y,z)) dV 
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Theorem: — JadV = J — dV + J v. (a v)dV 

Clt y y y 

Conservation of mass: a = p mass density 



AjpdV = = J^dV + Jv.(pv)dV 

UL v v u v 

| + v.(pv) = o = | + (v.v)p + P (v.v) = o 




Energy Stored (W E ) Rate Power Force Density 

Dissipated P E 



Work Rate = Mechanical Power 

- 1 1 - 1 2 

F = p f E — E Vs (force per unit volume) 

nt/m 3 
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f = fFdV 

t 1 

force (nts) 



D. Magnetoquasistatics 



V = J f , E' = E + vxB => J f ' = J f = ctE' = a(E + vxB) 
B = n(x,y,z)H 

jv.(ExR)dV = ^ExH-da" = -£V k I k = -jH.^-(n(x,y,z)H)dV 



-IP" 

v 

P dissipated =|E'.J' f dV = |E'.J f dV 

v v 

J f .(vxB) = -Jf .(Bxv) = -(jf xB). V 

H-— ji(x,y,z)H \ = -h ^ — Ui x, y,z H 

3t L v ' J |i(x, y,z) 3t L v ' 



vxB Of dV 



1 5 



H (x, y, z) 5t 



i^i 2 (x,y,z)|H| : 



5_ 
5t 



d_ 
5t 



|ji/(x,y,z)|Hr 



^J>rY7zT 



— li (x, y, z) H — 
2 I I 5t 



n(x,y,z) 



-n(x,y,z)|H|' 



1 ^VS^IH 2 9 

J'2 ^^x^zj at 



[n(x, y,z)] 



Aj n . dv = o^|i + (v.v)n = o (v-v = o) 



H.|^(x,y,z)H]= c ' 



n(x, y,z) H ' 



1 i-i 2 
--H Vu-V 
2\ I p 
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at 



-n(x,y,z)|H|" 



dV + P d issipated 



Energy density W^ 



V 



J, xb--|h| 2 vh 



dV 



Fm = force density 



Mechanical Power 



W M = j I n (x, y, z) |H| 2 d V , P djssipated = J E ' • Jf dV =j" E ' • Jf ' dV = j" o |e j 2 d V 



v ■ 



Total Magnetic Energy 



F M = J f xB-^|h| 2 V|i 



force density 



V. Compressible Media 

A. Electroquasistatics (EQS) 

Ohmic media: J' = aE' 

Polarization dependent on mass density (p) alone, electrically linear 

D = s(p)E 

EQS Galilean Transformation: 3 = oE + p f v 
Jv.(ExR)dV = $ExH.di = -£V k I k =-jE.| r [e(p)E] dV 



5t 



-jE.(aE + p f v) dV 

v 



1 5 1 2/ x|ef 

—rr E (p) E 

s(p)5tL2 ^1 I 
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5 


"i (p)|e| 2 1 


£2 (p)| E f 3 


( 1 1 


dt 


2 s(p) 1 1 J 


2 St 


vS(p)J 



i-|[e( P )i 



E = t- 



— 12 

E 



+1 3s( P ) 



5t 



5_ 
5t 



■e(p) E 



1 i|2 ge( P ) 

2 I I 5t 



5s(p) = as(p)5p * +v .(-U (Conservation of mass) 

5t 5p 5t 5t v F / 



de(p) _ de(p) 
fit 5p 



(-V.(pv)) 



-Zv k i k = -( 



at 



fe(p)|E| 2 



r 1 l-l 2 5s(p) 

dV- J - E — ^ dV 



5t 



f l|,2 5s(p) 
,2! 1 5t 



v 2 ' ' dp 



-Jv. 



13 8 |=|2 - 

2 5p 1 1 



J a |e| 2 dV - J Pf E-v dV 

v v 

(pv)dV 

dV + Jpv^V 



-IeI 2 — 

2 I I dp 



dV 



1 5 s i-|2 , r — I 

— p — E v»n da+ v- v 

2 op I I J I 



1 5 8 |-|2 

2 p 5p-rl 



i|E| 2 ^V P dV 
2l I 5p J 



Iv k i k = J 



at 



-e(p)|E|" 



dV + J a |e| 2 dV 



1 5 s i-|2 , 

-p — E v«n da 

s 2 P 5pl I 



+ fv. 

v 



Pf E-||E| 2 Vs + V 



1 5 s i=|2 
— p E 

2 5p I I 



dV 
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where 



58 

dp 



Vp = V8 



, electric energy 



W E = J ±s (p) |E| 2 dV , P dissipated = j a |E| 2 dV (power dissipated) 



F E - Pf E-||E| 2 Vs + V 



1 3 s i-|2 

2 p 5p-rl 



force density 



r 1 3 s i-|2 , 

cb — p — E v«n da = 
J s 2 3p 



because as S -> », E da -> 



3W • - - 

ZV k I k =— ^ + P dissipated +{F E .v dV 



dt 



v v ' 

Mechanical Power 



B. Magnetoquasistatics (MQS) 

MQS Galilean Transformation: Jf ' = Jf, E' = E + vxB 

B = ji(p)H 

jv.(ExH)dV = cfExH.da~ = -XV k I k = _JR.|-[n(p)H] dV 



Pdissipated = J E'.j; dV=fE'.Jf dV 

V V 

Jf • ( V x B) = -Jf • (B x v) = - ( J f x B) • V 



-J[E'- vxBj.Jf dV 



H.|^(P)H] = 



H(p) dt^W' 



1 d 



|i(p) dt 



^ 2 (p)|Hf 



d_ 
dt 



1 ^ 2 (p) |h| 2 

2 n(p) I I 



^ 2 (P)|H 



2 e_ 

dt 
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d_ 

at 



-h(p)|h 



2 I I dt 



_d_ 
dt 



jn( P )dv = o => M£i + v.[n(p) v] = o 

Mp) = Mp)3p ■ ^ + V .(pv)=0 

5t 5p dt dt v ^ 

"(-V.(pv)) 



dt dp 



-Iv k i k =-J 



at 



^(p)|h| 2 



dV -|i|H| 2 ^MdV P H 



dt 



-J(3f xBj.v dV 



v 2 ' ' dt {2\ I dp V > 



-Jv. 



1 3|i i— 12 - 

— H pv 

2 3p I I 



dV + Jpv-V 

v 



-|Rl 2 ^ 



dV 



— p — In I v • n da + f v • <! V 



1 n ^ iHl 2 

— p — H 



2 I I 



IX W 



at 



-h(p)|h 



dV + P diS5 -^p^|H| 2 v.n da 

s 5 P 1 1 



+j" V« 

V 



- - 1 i-i 2 11 du i-i 2 

Jf xB H V|i + V -p— H 

211 v 2 dp'' 



dV 
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where 

9|i 



p Vp ~ V|i ' ma 9 netic ener gy 

W M = {^(p)|h| 2 dV, P dissipated = jE'.i dV =|E'.J f dV Power dissipated 

F M = J f x B - - |R| 2 V|i + V [ - p ^ |R| 2 1 force density 



f 1 Dli i — v- — — I — y- 

d) — p — H v«nda = because as S ^ oo, H da -> 

I 2 dp i i II 



211 U 3p 



-|2 — — I — 12 



= ir + p ^ +|fm.v dv 

Mechanical Power 



C. Conclusions 
Force densities 



EQS: F E = Pf E-||e| 2 Vs + V 



1 3s i-|2 
— p — E 



MQS: Fm = J f xB-||h| 2 Vu. + V 



ip^lEl 2 
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